Exclusive Higgs boson production at the LHC: 
hard rescattering corrections 



J. Bartels a) : S. Bondarenko a) \ K. Kutak a ) ,b) * and L, Motyka c )'< 

a ) // Institute for Theoretical Physics, University of Hamburg, Germany 
b ' H. Niewodniczanski Institute of Nuclear Physics, Krakow, Poland 
c ) DESY Theory Group, Hamburg, Germany 
d ) Institute of Physics, Jagellonian University, Krakow, Poland 

January 16, 2006 



Abstract 

We examine, as a correction to the central exclusive Higgs boson production in pp collisions 
at the LHC, the rescattering of gluonic ladders off the proton. As usual, at the lowest order the 
hard part of this process can be described as a fusion of two hard gluonic ladders. We calculate 
corrections to this hard amplitude which are due to rescattering of these ladders. These corrections, 
which contain high mass diffractive excitations of the proton, have not yet been taken into account 
by the usual soft survival probabilities. Wc find that the correction due to the exchange of a single 
hard rescattering is negative, large and infrared sensitive. As a first step towards a more reliable 
description we therefore replace the rescattering exchange by a unitarized amplitude using the BK 
equation which generates the saturation scale Q s {x). We also include a soft gap survival probability 
factor. We discuss the results and outline possible future strategies. 
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1 Introduction 



One of the main experimental goals to be achieved at the Large Hadron Collider (LHC) is the discovery 
of the Higgs boson and the investigation of its properties. The discovery potential of LHC should be 
sufficient to find the Standard Model Higgs boson as well as the supersymmetric Higgs particles in the 
wide range of masses, starting from the lower limit of 114.4 GeV imposed by the LEP data up to 
the mass of a few hundred GeV. Precise electroweak measurements indicate, however, that the mass 
of the Higgs boson should be rather at the lower end of the allowed window, with the central value at 
113 GeV being determined from radiative corrections, and the upper limit of about 240 GeV at the 
confidence level of 95%. It has been proposed that, if the Higgs boson mass, Mh, is close to 

the expected central value, a measurement of the exclusive diffractive reaction should be possible: 



with the Higgs boson separated by rapidity gaps from the protons. The main advantage of this process 
is an extremely low background which is much smaller than the potential signal. The background 
suppression permits precise measurements - for instance of the quantum numbers of the Higgs boson. 
Besides that, in the exclusive production process, detection of protons provides a complete information 
on the invariant mass of the produced system. This opens an exciting possibility to find resonances 
corresponding to new heavy particles that would be invisible in conventional measurements, see |161 
1171 \1H\ . The main drawback of the exclusive Higgs production channel, however, is its very low rate. 
The requirement of leaving the protons intact and of producing no secondary particles reduces the 
cross section by many orders of magnitude. In the series of papers 0-^5 the authors discussed 
extensively the theoretical and phenomenological aspects of the exclusive production of heavy particle 
in pp collisions in perturbative QCD, and their approach will be used as the basis of our study. 

The mechanism of the considered process is illustrated in the diagram in in Fig. ^i. Thus, the 
hard subprocess is given by a fusion of two colour singlet gluonic ladders into the Higgs scalar. The 
diagram may be evaluated in the kt -factorisation framework using the unintegrated gluon densities. 
The required absence of hard radiation leaves un-compensated the virtual corrections to the process. 
Those corrections turn out to be large as they contain the scale of the Higgs boson mass generating 
large double and single logarithmic contributions. The necessary resummation of these logarithmically 
enhanced terms is realized by including the Sudakov form factor, which not only reduces the cross sec- 
tion but also provides the infrared stability of the hard exclusive scattering amplitude. Unfortunately, 
the basic hard subprocess is accompanied by multiple soft rescattering, which, typically, is inelastic 
and suppresses the exclusive cross section. Currently, the soft inelastic rescattering may be estimated 
only in phenomenological models |191 1201 1211 122j . Present calculations of the gap survival probability 
rely predominantly on a two channel eikonal model for the pp opacity, incorporating effects of the 
low mass diffraction and of the pionic cloud. The model was developed in Refs. [201 121j . and the 
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Figure 1: The exclusive Higgs boson production in a pp collision: a) the Two Pomeron Fusion contri- 
bution and b) a hard absorptive correction. 

parameters were fitted to describe the data on the total pp and pp cross section and on the elastic 
pp and pp scattering. 

In this paper we shall address the problem of computing a hard absorptive correction to the 
exclusive production process, which has already been mentioned in the literature JH] but has not been 
studied in detail yet. In brief, this contribution comes from an absorption of the hard gluonic ladders 
that fuse to produce the Higgs boson in the pp scattering process, see Fig. This contribution can 
also be viewed as the rescattering of a large mass intermediate diffractive state moving along one of 
the protons (the proton at the top of in Fig. ^>; let us call it the projectile) off the other proton (the 
proton at the bottom of in Fig.^3, the target), and this contribution is not part of the soft rescattering 
of the spectators that contributes to the standard opacity. In Fig. ^5, the large mass diffractive state is 
indicated by a dashed line. In perturbative QCD, within the leading logarithmic 1/x approximation, 
the rescattering is driven by the exchange of a BFKL Pomeron |23| I24 | 1231 126 which couples to the hard 
gluonic ladder by the triple Pomeron vertex j27J |2H1 HH1 EDI • Such a contribution is potentially large, 
since the BFKL rapidity evolution between the target proton and the triple Pomeron vertex extends 
over the length of Y > 14 and thus may generate an enhancement by two orders of magnitude. If this 
correction is really as large as anticipated, it is clear that also higher order unitarity corrections to 
this hard rescattering will be sizable and should be evaluated as well. 

As we shall discuss in detail, the situation is even more complicated. Namely, when evaluating the 
rescattering correction based on a simple BFKL Pomeron exchange, it turns out that the amplitude 
- which is infrared finite - exhibits a high sensitivity to the details of the infrared region. To be 
more precise, the momentum scale around the triple Pomeron vertex tends to be small. At first sight, 
therefore, this naive approach seems to be highly questionable. However, if we include - perturbative - 
unitarity corrections, the problems with the infrared sensitivity may be greatly reduced. As a guiding 
example let as quote results of Refs. [341 135j on the impact of a special class of unitarity corrections on 
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the behaviour of the gluon distribution at low momenta. In these references the Balitsky-Kovchegov 
(BK) equation [31 ( 132 1 133] was investigated, as a model which unitarizes the BFKL Pomeron exchange 
by a resummation of the BFKL Pomeron fan diagrams. It was found that the equation generates a 
momentum scale, the saturation scale Q S (Y), which provides an effective lower cut-off on the gluon 
momenta. The saturation scale grows exponentially with the evolution length Y, Q S (Y) ~ exp(AY), 
with A ~ 0.3, and the saturation scale enters the region in which perturbative QCD is applicable if Y 
is sufficiently large. Thus, in the presence of the perturbative saturation scale the high sensitivity 
to the infrared domain is eliminated. Let us stress here that the concept of a rapidity dependent 
saturation scale is rather robust. The concept has received strong support from phenomenological 
considerations, being the central part of the very fruitful saturation model [36[ 137] . It was also proven 
that the generation of the saturation scale is a universal phenomenon in a wide class of non-linear 
evolution equations [SHI OH HHj , which gives some confidence that the emergence of the saturation 
scale is not an accidental feature of the BK equation. 

In the second part of our investigation, therefore, we shall replace the linear BFKL equation for 
the screening Pomeron by the non-linear Balitsky-Kovchegov evolution equation. This corresponds 
to the resummation of the BFKL Pomeron fan diagrams originating from the target proton: this 
resummation will be shown to stabilize the correction amplitude in the infrared region. Certainly, in 
doing so we have made rather crude approximations. Namely, we have not taken into account other 
important unitarity corrections present in the case of proton-proton scattering, for instance those 
related to possible interactions of the Pomerons in the fan diagrams with the projectile proton (for 
a more detailed discussion see Section 7). For the total cross section of nucleus- nucleus collisions 
an effective field theory of interacting BFKL Pomerons has been formulated [411 1421 143j . For pp 
scattering an analogous framework is still to be developed. We therefore view our study only as a first 
step towards constructing a theoretical framework for evaluating hard absorptive corrections to the 
exclusive high mass production in pp collisions. We shall compute the magnitude of the correction in 
this, rather crude, approximation in order to provide some insight into how important the discussed 
effects might be. 

The structure of the paper is the following. In Section 2 we briefly review the standard framework 
introduced in the series of papers [5]-|llj. We derive the hard rescattering correction in Section 3. The 
calculation of the soft gap survival probability for the hard corrected amplitude is given in Section 4. 
In Section 5 the Balitsky-Kovchegov evolution of the gluon density is described. The results of the 
numerical evaluation of the hard rescattering corrections and the obtained cross sections are presented 
in Section 6. Implications of the results and future strategies are discussed in Section 7. Finally, 
conclusions are presented in Section 8. 
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2 The hard amplitude for exclusive Higgs production 

In this Section we briefly review the main features of the standard QCD approach to the hard 

part of the exclusive Higgs boson production in pp collisions. 

The basic diagram is given in Fig. the leading contribution of this amplitude is imaginary. 
The diagram is evaluated in the high energy limit, where the t-channel gluons are in the so-called 
nonsense polarisation states, and their virtualities are dominated by the transverse components of 
their momenta. The matrix element describing the fusion of two gluons into the Higgs boson is given 
by the top quark triangle, taken in the limit of a point-like coupling. This effective approximation 
is very accurate as the gluon virtualities are much smaller than the top quark mass. The amplitude 
of finding, in the proton, two gluons that carry tranverse momenta k and k and fractions x and x' 
of the longitudinal momentum of the parent proton is represented by the off-diagonal unintegrated 
gluon density f° (x,x';k,k'; fi), assuming that the density is probed at the scale \i. The three sub- 
amplitudes (Fig. 12), corresponding to the gluon emissions and to the Higgs boson production, are 
convoluted using the /^-factorisation. 

The kinematics of the exclusive Higgs boson production determines the arguments of the off- 
diagonal gluon densities. Thus, the elastic form-factor of the proton provides a sharp cut-off, which 
we assume to be Gaussian: ~ exp[— (k — k') 2 R 2 /4] with R 2 ~ 8 GeV 2 , on the momentum transfer 
carried by the gluonic system. Therefore, it is appropriate to write 



For future use, let us also define the unintegrated gluon distribution depending on the position b in 
the transverse plane 



f°^(x, x ; k, k ; /i) ~ f°^(x, x ; k 2 ; 



H) exp[-{k-k') 2 R 2 /4}. 



(2) 



ff(x,x , ;k 2 ,b-^) = J 




e^ b ff(x,x';k,q 



k;fj,) 



(3) 



(27T) 2 
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where q = k — k'. 

The next point that should be discussed is the disparity between the longitudinal momentum 
fractions x and x' . The production of the heavy scalar particle occurs close to the central region 



where x ~ Mh/^/s and x' ~ y/k 2 /s, so i' C i. At the leading logarithmic level the off-diagonal 
(w.r.t. x) distribution reduces to the diagonal one. Beyond the LL accuracy the effect of different x 
values of the gluons may be approximated by a constant factor ^IJ 

2 2 ^ r(A + 5/2) 

^ = r(A + 4) ' (4) 

provided the diagonal gluon distribution has a power-like x-dependence, f g (x,k 2 ; n) ~ For 
moderate scales one has A ~ 0.25. 

The last important point is the scale dependence of the unintegrated gluon density. In the BFKL 
approach to the evolution of the gluon density there is no dependence of the gluon density on the 
external scale. This happens because the virtual perturbative QCD corrections to the coupling of 
the ladder to the Higgs production vertex generate only the DGLAP logarithms log(/^ 2 /A; 2 ) (double 
and single) and not the small x logarithms. In double logarithmic (in log(/i 2 / k 2 )) accuracy, the large 
corrections may be resumed into the Sudakov form-factor 

V Jk 2 V 7T Jp UJ J 

The calculation of the Sudakov form- factor at single logarithmic accuracy has also been performed , 
giving 



T g (k,fi) = exp ( - / dzz 



(6) 



where 5 = is chosen to provide the correct angular ordering of the real gluon emissions, and P ( 
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and P qg are the DGLAP leading order splitting functions. For the exclusive Higgs production the 
scale /i ~ 0.62M# should be chosen in order to match exactly the first order single logarithmic 
correction. In our calculation we shall adopt the scale [i = M/^/2, which does not introduce a 
significant change of the result compared to the choice of . 

A thorough discussion of how to determine the two-scale diagonal unintegrated gluon distribution 
has been presented in Refs. |45| I4fi| . A simple and economic way to incorporate the Sudakov form- 
factor into the unintegrated gluon distribution is given by the approach in which the last step of 
the DGLAP evolution is resolved in terms of the transverse momentum. In this approach, the two- 
scale unintegrated gluon distribution is obtained from the collinear gluon distribution in the following 
way EH Enillll: 



f g (x,k 2 ;fi) = Q 2 



xg(x,Q 2 ) ■ T g (Q,n) 



(7) 

l 2 =k 2 



dQ 2 

In the off-diagonal case which is relevant for our present study, the logarithmic virtual corrections are 
associated with the production vertex at which the pair of gluons merges (see Fig. [^J), in contrast to 
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the standard case of a diagonal gluon distribution in DIS or in a hard inclusive production processes. 
Thus, in the Sudakov form-factor has to be replaced by its square root. 



Collecting all these results we arrive at: 

ff (x,k 2 ^) = R ( Q 2 ,J 



xg(x,Q ) • JT g (Q,n) 



(8) 

Q2 =fc 2 



dQ 

After contracting the polarisation tensors of the gluons with the ggH vertex we obtain the dominant 
imaginary part of the hard exclusive amplitude for pp — > pHp in the forward direction in the 
following form: 

ImM (y) = 2vr 3 A J ^ ff (xi, k 2 ; ff(x 2 , k 2 ; M ). (9) 



Here the constant A is given by 

A 2 =K V2G F a 2 s {M^) /9 vr 2 , (10) 

where Gf is the Fermi coupling, and K ~ 1.5 is the NLO i^T-factor. The Higgs boson rapidity in 
the centre of mass frame, y, defines the values of longitudinal momentum fractions: 

xi,2 = —7^ exp(±y). (11) 

The differential hard cross-section following from this amplitude takes the form: 



dy dt\ dt 2 



t „ b ,o 256^3 



After the inclusion of the proton elastic form-factors in the non-forward direction (parameterized by 
~ exp[— (k — k') 2 R 2 /4\ in @) one obtains the dependence of the cross section on t\ and t 2 of the 
form of exp(i? 2 (ii +i2)/2). Thus, after the integration over momentum transfers t\ and the single 
differential cross section reads: 

da P p->pHp(y) _ \M (y)\ 2 

dy 64ir 3 R 4 ' 1 ; 

Note that, in Q, the form of the integral over the transverse momentum suggests an important 

sensitivity to the infrared domain. The rescue comes from the Sudakov form- factor which suppresses 

the low momenta of the unintegrated gluons faster than any power of k. It turns out that the dominant 

contribution to the integral comes from momenta k > 1 GeV. 

More details about the derivation of the amplitude of the exclusive Higgs boson production are 
given in the series of papers |5]-|11|. In what follows we shall refer to the amplitude @ as to the Two 
Pomeron Fusion (TPF) Amplitude. 
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Figure 3: The triple Pomeron vertex. 
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3 Hard rescattering corrections 
3.1 Construction of the amplitude 

The amplitude © includes resummed QCD corrections up to single logarithms of the Higgs boson 
mass. The large available rapidity gap between the scattered protons (itot ~ 19 f° r the LHC energies), 
however, suggests that significant QCD corrections enhanced by small x logarithms should be also 
considered. The spectator rescattering is expected to be of non-perturbative nature, and commonly 
it is accounted for by a soft gap survival probability factor. 

Fig. |1] illustrates a process of absorption of the perturbative gluonic ladder from the projectile 
(upper) proton in the target (lower) proton by the exchange of a QCD Pomeron across a large rapidity 
gap. Let us stress that the evolution length in rapidity of the screening Pomeron is large only if 
it couples to the ladder between the projectile proton and the produced Higgs boson. To be more 
specific, for the central Higgs boson production at the LHC the rapidity distance between the hard 
ladder that originates from one proton, and the other proton, y^ = log(l/x4) > 14, is large while 
the rapidity range occupied by each of the two hard ladders, is much smaller, y\ = log(l/xi) < 5 
and 2/3 = log(l/x3) < 5. In contrast to this, the rapidity distance y^ can be much larger, and a 
small x (BFKL |2Ej) resummation of the perturbative QCD series may lead to a large enhancement 
of the amplitude ~ exp(Ai/4) with A ~ 0.3. Therefore this correction may be sizable, despite the 
fact it is suppressed by the factor of a 2 s at the triple Pomeron vertex, illustrated in Fig. El Thus, we 
shall evaluate the absorptive correction corresponding to the rescattering of the hard projectile gluonic 
system off the target proton. The coupling of the screening Pomeron to the screened gluonic ladder is 
known in the high energy limit; it is derived from the triple Pomeron vertex |27[ I28 | 129 } I30j . 

The triple Pomeron vertex in pQCD has been computed in the LLl/x approximation |27l I28j . 
and it follows from the 2^4 gluon vertex by projecting on the colour singlet states of the diagram 
represented symbolically in Fig. |3J In the process of interest the momentum transfers of all three 
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Figure 4: The hard rescattering correction to the exclusive Higgs boson production. 



Pomerons are small: for the upper Pomeron and for the lower left one they are damped by the 
elastic proton form- factors. We, therefore, will neglect them in the expression for the triple Pomeron 
vertex, and it is will be sufficient to evaluate the vertex in the triple-forward limit. In this limit, the 
convolution [V3p<8>-D] of the triple Pomeron vertex V^p with the two-gluon amplitude above the vertex 
(here we follow the notation of j28j . below we will replace D by the unintegrated gluon density) is 



where ±&2 an d ±fc 3 are gluon momenta inside the two Pomerons below the vertex (see Fig. OJ). Note 
that, in formula (|14jl. the triple Pomeron vertex includes both real (the first term in the l.h.s.) and 
virtual (the next term in the l.h.s.) contributions. The normalisation factor, c 3 p, of the vertex will 
fixed by the form of the final expression for the correction amplitude (|18[) . and we have constructed 
this expression in such a way that c 3 p = 1. 

In order to build up the complete expression of the amplitude of the diagram shown in Fig. 0] it is 
necessary to take into account the following building blocks: 

1. The unintegrated gluon distribution, fg {x\, k\ ; /i a ), describing the amplitude of the gluons 



given by [23 EH]: 




(fci + fc 3 ) 2 1 (fci + fc 3 ) 2 + k\~ (fcx - k 2 y + (fcj + fc 3 ) 2 



(14) 
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above the triple Pomeron vertex V%p. This gluon distribution is evaluated at relatively large 
x\, and the momentum scale in the vertex should be given by the typical virtualities of the 
vertex, fj, a ~ k%. So the gluon distribution will be evaluated using Q , and we shall denote the 
convolution (|T1)) of the triple Pomeron vertex with the unintegrated gluon, fg(x\,k\\k\), by 

[VsP^f^HxxM^h). (15) 

2. The two-gluon Green's function Ga(&2, k 3 ; x±, x 2 ; £i , [i b ) that describes the propagation of the 
two gluon singlet between the Higgs production vertex and the triple Pomeron vertex. Two ap- 
proximations for the Green's functions will be considered: the two gluon approximation {Gq{. . .)) 
and the diffusion approximation to the BFKL evolution (G*bfkl(- • •))• The Green's function will 
have the typical scales fi a and /if,, where the latter scale belongs to the Higgs production vertex. 
As to the scale of the triple Pomeron vertex, [i a ~ A>2, a priori, we do not know whether it will 
be large. However, the large scale is certainly present in the Higgs vertex, [x b ~ Mh/2. In order 



to account for the scale dependence at the Higgs vertex, the Sudakov form factor, y T g (k 3 , fj, b ), 
will be included into the Green's function: 



G (k 2 ,k 3 ;x 1 ,x 2 ;Hb) = yT g (h,[i b ) k^5 {2) (k 2 - k 3 ) (16) 

and 

i°g 2 (fc 2 2 A§) 



n „ , s k 2 k 3 ^T g (k 3 ,fi b ) 

G B FKL{k 2 ,k 3 ;xi,x 2 ;n b )= — Y = exp 

2vr ^7ri;iog(xi/x2) 



AD log(xi/x 2 ) 



, (17) 



where only the leading conformal spin was taken into account, ujq = 41og(2)a s is the LL BFKL 
intercept, D = I4((3)a s is the diffusion coefficient, and a s = 3a s /n. 

3. The Higgs boson production vertex which does not differ from the TPF case. 

4. The off-diagonal unintegrated gluon distribution, f g b \x^, k^,q — k 3 ; fi b ), of the gluon distribution 
that enters the Higgs boson production vertex from the lower side, with the total transverse 
momentum q. The gluon distribution is taken to be the same as in the case of the Two Pomeron 
Fusion, i.e. it is an off-diagonal gluon distribution given by formulae (JEJ and (J2J) and with 
(i b = M H /2. 

5. The off-diagonal unintegrated gluon distribution, f g c \x^,k^,—q — k^;fi c ), for the screening 
Pomeron that propagates across the large rapidity distance from the lower proton to the triple 
Pomeron vertex with momentum transfer — q. The scale [i c is set by the gluon virtuality, so 
one has f g (£4, k^, —q — k^) = f g c \x^, k^, —q — k^; k^). For a detailed discussion of this gluon 
distribution see Section 4. 



1 Note that, for fi a — k\, the Sudakov form factor, T g (ki, fi a ), in Q introduces only a minor effect. 
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After these definitions the expression for the dominant imaginary part of the diagram of Fig. 0] 
can be written in the following form: 

rx b dx 4 f d 2 q f d 2 k 2 f d 2 k 3 f d 2 k 4 



„,m , n 9 n . , f x b dx 4 f dr q f d z k 2 f d z k 3 f 



2ixk\ 



x [V3P®/i a) (^i)](fc2, fc4)GA(fc2,fc3;xi,x 2 ;M H /2) 



x f^(x 3 ,k 3 ,q-k 3 ;M H /2)f^(x 4 ,k 4 ,-q-k 4 ), (18) 

where | = tends to unity in the large N c limit, and the coefficient 2tt 3 A has been extracted in 

order to enable an easier comparison with the Two Pomeron Fusion contribution given by ©. x a and 
xj, denote the lower and the upper limits of the integration over x 4 , respectively. The fractions of the 
longitudinal momenta of the upper proton x\ and x 2 and of the lower proton x 3 , x 4 read: 

M H M H k 2 

x 3 = exp(-y), x 2 = exp(y), x\ ~ (19) 

ys v s X 4 S 

k 2 k 2 

x a <x 4 <x b , £ a ~-^, x b ~ — . (20) 

s X 2 S 

The energy scale feg that enters the definitions of gluon xi corresponds to typical values of the vir- 
tualities, and we set Uq = 1 GeV 2 . The dependence on the momentum transfer g, circulating in the 
Pomeron-proton triangle, appears only inside the off-shell unintegrated gluon distributions, fg^ and 

(c) 

fg , which contain the rapidly varying proton form-factor. The q dependence of the other vertices 
and of the Green's function, Ga, is weak and will be neglected. The amplitude (|18ft includes a factor 
of 2 coming from the possible emissions of the Higgs boson from both of the Pomerons below the triple 
Pomeron vertex. In addition, we have to add another contribution in which we interchange the role 
of projectile and target: since the structure of both protons is identical, we simply replace y — > —y. So 
finally the correction amplitude may be written as: 

M COII (y) = M« r (y) + M^ T (-y). (21) 

In the case of the central production, y = 0, and M corr (y) = 2Mco} r (y). 



3.2 Infrared stability 

It is instructive to study the infrared behaviour of the integral in (|18|) using a simplified framework and 
the saddle point method. First of all, we note that the integrals in (|18j) are infrared finite. This follows 
from properties of the triple Pomeron vertex (which vanishes whenever one of the external momenta 
goes to zero) and of the unintegrated gluon densities (which, as a result of gauge invariance, also 
vanish as one of the two gluon momenta goes to zero). Next, after averaging over the angles in k 2 and 
k 4l the triple Pomeron vertex (|T1|) imposes the conditions k 2 > k\ and k 4 > k\. Finally, disregarding 
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irrelevant constants and details one finds that, for the real contribution of the triple Pomeron vertex, 
the correction term (|18jl behaves as 



where we used, for the Green's function Ga, the approximation of the two-gluon exchange, see (|16|) . 
Note that the Sudakov form-factor, \/T s (k2 , Mjj/2), is present in the Green's function (shown explic- 
itly in l)22[l) and, implicitly, in the gluon distributions fg b \xs, k^) Mh/2). 

In this analysis we have introduced a cut-off, k$, of the low momenta in the integration over k%. 
The scale ko is meant to lie in the perturbative domain above the low hadronic scale of the proton 
form- factor. The structure of the nested integrations in (|22|) implies that this scale ko is a lower limit 
for all momenta. By studying the variation of the integrals (|22j) as a function of ko, we obtain a 
measure of the sensitivity of the correction to poorly known details of the proton structure at low 
momentum scales. The strong coupling constant, a s , is kept fixed, as the running coupling introduces 
only logarithmic corrections, and, for the moment, we are interested in tracing only the leading power 
of ko that arises from the integral 1)22(1 . In our estimate we use the Sudakov form- factor in the double 
logarithmic approximation: 



Let us start by performing the analysis of Qr with the conservative assumption of a vanishing 
anomalous dimension of the gluon density (i.e. a constant behaviour of the gluon density as a function 
of momentum). Then the integral over &2 behaves in the same way as the integral defining the Two 
Pomeron Fusion amplitude. It is dominated by a contribution from a saddle point at a perturbative 
value, k2 = (Mh/2) exp(— 2ir / 3a s ), with only marginal sensitivity to the infrared domain. 

The nested integrations over k\ and k^ in (|22j) are dominated by momenta close to the cut-off 
scale ko and lead to a result ~ l/^o- This means that our correction is not small, but the magnitude 
is sensitive to details of the gluon distributions f g a ^ and fjf^ at low momenta, and hence, unreliable. 
At first sight, this seems to imply that the whole Ansatz of choosing, for the absorptive correction, a 
perturbative gluon ladder, cannot be justified. 

A closer look, however, shows that this may not be the final answer. Namely, many arguments 
have been given that the gluon density at very small x should saturate, i.e. it scales and depends upon 
the ratio k 2 /Q^(x). Here Q s (x) ~ x~ A is the saturation scale which grows in 1/x with A ~ 0.3. This 
leads to a substantial change of the A; 2 -dependence of the gluon density for momenta being smaller 
than the saturation scale: the saturation scale becomes an effective lower cut-off of the gluon momenta 
in the small-A; 2 region. The phenomenological analysis of HERA data |361 I37| indicates that Q 2 (x4) 
takes values of a few GeV 2 for the range of x± relevant in our problem where log (1/2:4) can be as large 
as 14 to 19 units. 






(23) 
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Applying these arguments to our integral (|22[) we find that, in the presence of the saturation scale, 

(c) 

the momentum dependence of gluon distribution fg is modified, and the nested integrals over k\ and 
&4 give 



With Q s {x^) being well above the infrared scale k$, this represents a well-behaved result with a low 
sensitivity to the infrared physics. Keeping in mind that, in ([22)1. the integration over ki is infrared- 
safe, one finds that Qr is infrared stable. 

An analogous reasoning can be applied also to the contribution of the virtual correction part of 
the triple Pomeron vertex to the integrals in (|18j) . The behavior is approximated by: 



Qv ~ J k2 ^ J k2 ^ log(k 2 2 /kl)f^(x 2 , k\) fW(x 3 , k\; M H /2) if(x 4 , k\) ^T s (h , M H /2), (25) 



and a similar term, Q' v , is obtained by interchanging f g (x^, k\] Mh /2) — ► /g (#3, fc|> Mg/2) and 
fg C \xi,k2) — > fg C \x^,k1). One finds that the integral over &2 (&i) in (Qv) is stabilized in the 
infrared by the saturation scale, and the integral over k\ (A^) is stabilized by the Sudakov form-factor. 

Collecting the estimates for Qr, Qy and Q' v one concludes that the perturbative determination 
of the amplitude of screening correction is possible provided that the gluon distribution fg C \xi,k\) 
possesses the saturation property, i.e. the simple gluon ladder that we have used initially has to be 

(c) 

corrected. In this case the momentum scale at the upper end of the screening correction, f g , is 
not small, and our perturbative Ansatz is justified. Note that, in this context, saturation plays a 
role quite analogous to the Sudakov form factor at the production vertex: without the Sudakov form 
factor the momentum integral near the production vertex would not have been infrared-safe. Similarly, 
saturation renders the nested k\ and k^ integrations infrared-safe. 

4 Nonlinear evolution of the gluon distribution 

In this section we use the simplest available tool, the nonlinear Balitsky- Kovchegov (BK) equation |31[ 
1321 133| . in order to model the unitarity corrections for the screening Pomeron ladder in Fig. |I] This 
equation resums fan diagrams of BFKL Pomerons (Fig. in the LLl/x approximation, and it has 
been derived, in the large iV c limit, for the case of a small probe scattering off a large nucleus. The BK 
equation generates, in a natural way, the saturation scale growing exponentially with rapidity |291l35j . 
The equation has originally been derived for the dipole scattering amplitude, but later on it was re- 
formulated as an explicit equation for the evolution of the unintegrated gluon distribution |481 149j , 
and in this paper we employ the latter form. 

The hard rescattering correction amplitude in (|18|) was given in the momentum representation, in 
order to match the variables of the triple Pomeron vertex. In this section, however, we find it more 




(24) 
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Figure 5: The fan diagram approximation to the hard rescattering correction to the exclusive Higgs 
boson production. 

convenient to do a partial Fourier transform from momentum transfer carried by the Pomerons to the 
impact parameter variable. The impact parameter space seems to be advantageous, both in performing 
the non- linear gluon density evolution and for evaluating of the soft gap survival factor. The correct 
treatment of the impact parameter dependence of the scattering amplitude is difficult in the BK 
framework. In the original formulation of the BK equation, the gluons are assumed to be massless, in 
accordance with perturbative QCD; this contradicts, however, essential features of non-perturbative 
QCD dynamics leading to the confinement of colour. In particular, in the perturbative framework 
one finds contributions from large dipoles which generate a power like tail of the impact parameter 
profiles |5U1I51| . Such a tail leads to a violation of the Froissart bound [52] • Therefore, the limit of large 
impact parameters of the BK must be modified, using some phenomenological model. In practice, all 
the applications of the BK equation to the proton structure are based on the approximation of treating 
the evolution at fixed impact parameter |35 | I48[ 153* 1 I54j . i.e. the evolution proceeds independently at 
each value of impact parameter b. 

Thus, our starting point is the form of BK equation proposed in Ref. |48| : 

df g (X,k 2 ) _ a S N Ck2 f dk' 2 I fg (X, k' 2 ) - fg (X, | fg ( X , jg) l 

<91ogl/x vr Jk'2 k' 2 \ \k' 2 -k 2 \ [Ak' A + k^ } 

2 

(26) 

where the assumption has been made that the distribution of gluons is uniform inside a transverse 
disc with the radius R. The unintegrated gluon density in the impact parameter plane therefore takes 
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>A\ 2 — 

dk 2 ) R 2 



at 1 - k 



00 dk' 2 k' 2 \ 
k2 l^ l0g [^ f ^ k 



/2\ 



the form: 

f g (x,k 2 ,b) = f g (x,k 2 ) 9 -^^. (27) 
In 1)26(1 . the factor 1/R 2 in front of the non-linear term is a result of the impact parameter integration: 

d 2 bf g (x,k 2 ,b)f g (x,k 2 ,b) = -L f g (x,k 2 )f g (x,k 2 ). (28) 

In our case, however, it is necessary to control the impact parameter dependence in more detail. 
Therefore we assume, at xq = 0.01, an initial Gaussian impact parameter profile: 

/,(*„, k 2 , b) = f g (x , k 2 ) eXp( ^ 2) , (29) 

and we set the radius R ~ 2.83 GeV -1 (R 2 = 8 GeV -2 ), in accordance with the proton radius seen in 
the i-dependence of the exclusive J /ip production at HERA. The evolution equation for the density 
of the unintegrated gluon reads: 

dfg(x,k 2 ,b) _ asNc^ f dk' 2 f fg (X, k' 2 , b) - fg (X, k 2 ,b) | fg (x , ^, b) \ 



<91ogl/x vr hi k' 2 \ \k' 2 -k 2 \ [Ak' 4 + k A }^ 



wji- k 2 4^) k 2 



2 



r^dk' 2 (k' 2 \ ~ a ' 



(30) 



roo dk' 2 , f k' 2 

dk 2 ) 

This evolution is local in b, but due to the non- linear term, for x < xq, the impact parameter 
dependence cannot be factorized out from the resulting unintegrated gluon density, f g (x,k 2 ,b). 

It is clear, however, that the BK equation (|26[) needs many improvements to be applicable in 
phenomenology, in particular it is essential to take into account, in some approximate manner, the 
NLLl/x corrections to the linear BFKL kernel |55| 156) . Following the path proposed in series of 
papers |57l 1481 H9"] . we include, in the linear evolution kernel, the following improvements: 

• the consistency constraint |58| I59j that models an essential part of the NNLl/x correction to the 
BFKL kernel and realizes, approximately, also a resummation of the corrections to all orders; 

• the non-singular part in the x — > limit of the DGLAP splitting function is added; 

• the DGLAP quark density evolution is also taken into account, and the evolution equation is 
coupled to the gluon density evolution using the DGLAP splitting functions; 

• the strong coupling constant is running with the scale given by gluon virtualities in the BFKL 
ladder. 

These modifications are consistent with the exact NLLl/x corrections to the BFKL kernel, and they 
impose the constraints on the evolution which follow from the DGLAP formalism in the collinear limit 
|60l 1611 16~2*] . The implementation of these corrections has been described in detail in the recent papers 



15 



|48l I49j . and the final form of the equation can be found there as well. It is rather straightforward to 
add, on top of the modifications described above, the impact parameter dependence, in a manner very 
similar to (|30|) . The normalisation of the input parameterisation for the gluon and quark distributions 
at high x has been tuned to provide a good fit of the F 2 data and the obtained distributions 
have been shown to give a good description of inclusive heavy quark production at HERA and at the 
Tevatron [51]. A detailed description of the properties of the obtained partonic densities will be given 
in a separate paper [63] . 

In our numerical calculations of the exclusive Higgs boson production l|18jl we will take f^ c \x,k 2 ,b) 
as being the solution to the collinearly improved, impact parameter dependent BK equation. Let us 
remind that the other two gluon distributions of p8jl. f^ a \x, k 2 , b) and f( b \x,k 2 ,b), are probed at 
rather large values of x ~ 0.01, and non-linear effects can safely be neglected. Instead, we apply the 
prescriptions and (jHJ) , using with the CTEQ5 parameterisation for the collinear gluon distribution 
functions. 



5 Soft gap survival probability 

Thus far, we have determined the hard component of the amplitude for exclusive Higgs boson produc- 
tion. The condition to leave the protons intact in the production process, however, imposes constraints 
on the soft rescattering of the protons. The soft gap survival probability factor, S 2 , quantifies the 
probability that the protons emerge from the reaction intact despite possible soft rescatterings. The 
gap survival probability can conveniently be estimated in the impact parameter representation by em- 
ploying, for instance, the two-channel eikonal model |20U21| . Therefore, we now translate the exclusive 
Higgs boson production cross section into the impact parameter representation. 

Let us start with a general discussion. We consider a scattering amplitude of an exclusive produc- 
tion in the impact parameter representation 

M(b l5 b 2 ) = lj^J j^M( qi ,q 2 )e^e^ b \ (31) 

where M(q 1 ,q 2 ) is the amplitude given as a function of the transverse momenta. Then, the integrated 
cross-section 

'-TZ/WU™*^ (32) 

can be expressed in the following way: 

a = ilb / dHl I d2b2 l^ bl ' b2 )l 2 ' ( 33 ) 

or, alternatively, 

° = lL / ^ I d% '* (61 ' b ~ bl)|2 ' (34) 
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Figure 6: The exclusive pp —>■ p H p amplitude with a hard and soft rescattering correction. 

where bi and &2 are transverse positions of the production vertex measured from centres of the protons, 
and b is the impact parameter of the collision. 

In the Two Pomeron Fusion case, the off-diagonal unintegrated gluon distributions are probed at 
moderate values of x ~ 0.01 and at virtualities k 2 > 1 GeV 2 , so one can assume the gluon density to 
take a factorized form: 

ff(x, x', k 2 , b; fi) = ff(x, x', k 2 ; M ) S(b) (35) 

with 

S(b) = ^- 2 e'^, (36) 

where R ~ 2.83 GeV -1 is the proton radius, measured in the exclusive production of J/ip production 
at HERA (cf. our discussion of the impact parameter dependent BK equation), and S(b) is normalized 



J d 2 bS(b) = 1 . (37) 



Therefore, the cross section for the production of the Higgs boson in the Two Pomeron Fusion ap- 
proximation takes the form: 

d^H P (y) = ±-Jd 2 bJ d'hlS^Sib - b 1 )M Q {y)\\ (38) 
with Mg(y) given by ©. Including the hard rescattering correction, we arrive at the following expres- 



sion: 



da^\y) pp ^ pHp = _1_ J d% J d 2 fei | 5(bi)5(6 _ 6i)MQ(2/) + i g corr(2/)5)6i) |2 ! (3g) 



dy 

where the dominant parts of both Mo(y) and M COTT (y, b, b±) are imaginary, and 

~ , , , . [9 o , o „. f x b dx 4 f d 2 k 2 f d 2 ks f d 2 k± 
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x [V3P® $ a) (xi,&-&i)](fc2, k A )G K {k 2 ,k- i ;x u x 2 -M H /2) 

xff\x^klbr,M H /2)ff{x^klb x )} - {y^-y}- (40) 

The impact parameter dependence of fg a ^ and fg^ is given by S(b) and may be factorized out, in 
analogy with the case of the Two Pomeron Fusion amplitude. The impact parameter dependence 
of fg is determined by the BK evolution, and it is different from S(b). Note that, in the 
convolution in the momentum transfer q leads to a product of fg b ^ and fg , evaluated at the same 
point, 61, in the transverse position space. 

Having defined M corr (y, b,b\), it is straightforward to write down the cross section for the exclusive 
production, taking into account the soft rescattering 



, (0+1), n , s. 



Jd 2 bJ d 2 h |5(6i)5(6 - 6i)M (y) + M corr (y, b, b^ 2 exp(-0( Sj b)), (41) 



dy 

where 0(s, b) is the opacity of the pp scattering at squared energy s and impact parameter b. In 
this paper we adopt the opacity factor of the two channel eikonal model proposed in [201 . The gap 
survival probability reads: 



pp — >pHp 

(y) 

dy 



S 2 = . (42) 

pp^pHp 1 

\ / y=0 

The impact parameter profile of the rescattering correction amplitude Q4UJI is different from the Two 
Pomeron Fusion case; thus the gap survival factor will take different values in those two cases. 



6 Results 

We have performed all numerical estimates of the amplitude and of the cross section for the central 
exclusive Higgs boson production, choosing y=0, assuming the LHC energy, yfs = 14 TeV, and taking, 
for the Higgs boson mass, Mjj = 120 GeV. As we have stated before, in this paper our main goal is 
a general understanding of the hard screening mechanism and an estimate of its magnitude, rather 
than a precise determination of the cross section. 

As a first step, we have considered the rescattering due to a single BFKL ladder. As expected, due 
to the long rapidity evolution of the BFKL amplitude this correction is large. As we have discussed 
before, since this single-ladder correction is highly sensitive to the infrared region, the precise numerical 
value of its magnitude is not reliable. From this part of our analysis, we therefore only conclude 
that this type of hard rescattering may potentially be large and, in any realistic estimate, unitarity 
corrections of the screening BFKL ladder have to be taken into account. 
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In the next step of our numerical analysis, therefore, we have made use of the BK-corrected 
unintegrated gluon density, fg. First, using (|4()jl. we have evaluated the hard rescattering amplitude 
as a function of the impact parameter, at fixed X4. The contribution of momenta lower than k$ = 1 GeV 
was assumed to vanish. The main reason for that is that the gluon densities are poorly known for 
gluon momenta lower than ko, so is the behaviour of gluon propagators, and in that region one has 
to rely on uncertain extrapolations. On the other hand, with the integrand of ()18|) being positive, 
the application of the cut-off leads to a conservative estimate of the correction term. Still, the main 
contribution to the correction amplitude comes from the range of momenta around the saturation 
scale Q s {xa)- 

Recall that, for the four-point gluon Green's function which connects the Higgs vertex with the 
triple Pomeron vertex, we consider two different scenarios - the simple two gluon exchange ()16|) and 
the BFKL Green's function in the diffusion approximation (|17|) . The strong coupling constant was 
fixed at a s = 0.15, both in the Green's function and in the triple Pomeron vertex. Such a choice is 
suitable for the LL BFKL Green's function, as the small a s reduces both the Pomeron intercept and 
the diffusion coefficient, in line with the NLO corrections to the BFKL evolution. This value of a s 
seems, however, to be too small for the triple Pomeron vertex: thus the estimate of the correction is 
conservative. Therefore, in the vertex, we shall consider also values of a s other than our default choice 
of a s = 0.15, while in the BFKL Green's function we always keep a s = 0.15. 

In Fig. [7|we present the total hard rescattering correction to the amplitude in the form of the ratio 

|(iM corr (y,b,bi)/dyi| 



Ry(yi 



(43) 

J b=b!=0,y=0 



\M Q {y,bM)\ 

with the impact parameter dependent Two Pomeron Fusion amplitude being given by (cf. Q and 

(El) 

M (y, b, bi) = M (y) S(b - b 1 )S(b 1 ). (44) 

Here Mo(y) has also been evaluated with the infrared cut-off ko = 1 GeV, and the off-diagonal 
unintegrated gluon distributions have been obtained from the CTEQ5 collinear gluon distribution. 
Note that we plotted the ratio as a function of the rapidity related variable, y\ = log(l/xi), instead 
of the less convenient variable X4 (the relation between these variables is given by (|19|1 ). The variable 
yi has the meaning of a rapidity distance of the triple Pomeron vertex from the projectile proton, and 
it is related to the invariant mass of the rescattering intermediate state by M\ ~ /cq exp(yi). 

In Fig. [3 we compare the results, obtained with the two choices of Green's function and with the 
two values a s = 0.15 and a s = 0.25. The relative sign of the correction term is negative. From this 
figure we draw the following conclusions: 

• The hard rescattering correction is the highest for the highest mass of the diffractive intermediate 
state. This is due to the strong suppression of the gluon density in the proton at large values 
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Figure 7: The ratio R y {yi) (see \4<ty ) °f the hard rescattering correction to the Two Pomeron Fusion 
term at b = 0, as a function of the rapidity distance, y\, of the Triple Pomeron Vertex from the 
projectile protons, for a s = 0.15 and a s = 0.25 at = 14 TeV. The black line belongs to 
the Green's function, Ga = GbfkL; the dashed line to the approximation of two elementary gluons, 
Ga = Go- 
of the gluon x. We interpret this as an indication that the contribution of the hard rescattering 
correction is clearly separated from the soft rescattering contributions of the two-channel eikonal 
model |20| l2T] which incorporates only the low mass diffractive intermediate states. 

• There is no significant difference between the BFKL Green's function and the exchange of two 
elementary gluons. The main reason is that the BFKL evolution length is short for larger values 
of yi, which give the dominant part of the corrections. The relative BFKL effects are more 
pronounced for small y\, but this region provides only a small contribution. 

• The relative magnitude of the correction is quite significant for a s = 0.15, and for a s = 0.25 the 
correction even exceeds the Two Pomeron Fusion amplitude. The interpretation of this striking 
result will be discussed in the next section. 



In Fig. |H1 results are shown for the impact parameter dependent correction amplitude, integrated 
over j/i, as a function of the transverse distance, b\, of the Higgs boson production vertex from the 
center of the target proton: 

\M COTT (y,b, b{) 



Rb(bi 



\M (y,b,h)\ 



(45) 



^=0,b=0 
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Figure 8: The ratio Rb(b\) (see j45\ )) of the hard rescattering correction to the Two Pomeron Fusion 
term at b = 0, as a function of the transverse distance, b\ of the Triple Pomeron Vertex from the 
centre of the target "proton, for the two values a s = 0.15 and a s = 0.25 at -^/s = 14 TeV. We only 
show the results for the BFKL Green's function. 

Here we have set 6 = 0, but, in fact, in the ratio the impact parameter dependence coming from 
the form factor of the projectile proton drops out, and any value of b can be taken. In Fig. |H] only 
the results obtained with the BFKL Green's function are displayed. There is not much difference 
between the BFKL Green's function and the two gluon approximation. One sees in the figure that 
the ratio of the absorptive correction profile, M COTT (y,b = 0,&i), to the Two Pomeron Fusion profile, 
Mo(y,6 = 0, bi), decreases with increasing distance, b\, from the proton center. Thus, the correction 
term tends to suppress the Higgs production inside the target proton. It is also clearly visible that the 
correction is large for both choices of a s . 

Performing, in (|41|). the integrations over the transverse positions, b\ and b, leads to the differential 
cross section of the exclusive Higgs boson production, daj^j^g (y) / dy which takes into account both 
the soft and the hard rescattering corrections. The results of this calculation, for various values of a s 
in the triple Pomeron coupling, are collected in Tab. ^ At a s ~ 0.2 the absolute value of the hard 
rescattering correction exceeds the Two Pomeron Fusion term, and this explains why, at this value of 
ct s , the corrected cross section takes its smallest value. For all values of a s that we have considered 
the cross section with the hard rescattering correction is smaller than the Two Pomeron Fusion cross 
section by factor larger than 2.5. Let us repeat that, by varying a s , we want to illustrate the sensitivity 
of the result to the size of the triple Pomeron coupling. Further discussion will be given in the next 
section. 
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da TPF /dy 


da /dy (a s = 0.15) 


da/dy (0.17) 


da/dy (0.2) 


da/dy (0.22) 


da/dy (0.25) 


0.4 


0.16 


0.12 


0.042 


0.08 


0.14 



Table 1: The differential cross section y apP dy PHp j Q fi n obtained from the Two Pomeron Fusion 
amplitude, and corrected by the hard rescattering and by the soft gap survival probability factor, for 
various values of a s at y/s = 14 TeV . 



S 2 (TPF) 


S 2 {a s = 0.15) 


S 2 (0.17) 


S 2 (0.2) 


S 2 (0.22) 


S 2 (0.25) 


0.024 


0.037 


0.061 


0.051 


0.035 


0.016 



Table 2: The soft gap survival probability factor, S 2 , for the Two Pomeron Fusion amplitude, corrected 
by the hard rescattering, for different values of a s at -y/s = 14 TeV . 

Numerical values of the survival probability, based on (|42[). are given in Tab. |2] for different values 
of a s . The variation of S 2 with a s is a result of the modification, due to the correction term, of the 
shape in impact parameter b of the hard amplitude. 

The definition of the survival probability, l)42j) . which we used is different from the definition of 
S 2 in |5]-|11|. |2()1 121j . In the case when we consider only the Two Pomeron Fusion term and the 
Gaussian form of the form factors, the definitions of S 2 in (|42jl and in Ref. [20] lead to the same 
numerical result. In the case of a general shape of the form-factors, however, we find it necessary to 
use the general formula for the soft gap survival factor of (|4U|1. (|4*T1) and (|42|). 

7 Discussion and future strategies 

The hard rescattering correction to the Two Pomeron Fusion amplitude, even after including uni- 
tarization of the single BFKL exchange, was found to be large. In particular, for a more realistic 
choice of the strong coupling constant, a s = 0.25, in the triple Pomeron vertex, the correction exceeds 
significantly the Two Pomeron Fusion term. As a consequence, we have to conclude that we have 



22 



found another potential source of a significant suppression of the exclusive Higgs boson production at 
the LHC. 

As we have said before, we cannot consider our results as representing a reliable numerical final 
answer. In our calculation we have considered only a sub-class of the screening corrections which - 
in a BFKL-based description of diffractive Higgs production - should be included, and in evaluating 
these diagrams, we have made severe approximations. In this paper, we have restricted ourselves to 
Pomeron fan diagram configurations in the exchanged screening system which lead to the Balitsky- 
Kovchegov evolution of the gluon density. Such a choice is a natural starting point for an estimate of the 
amplitude beyond the Two Pomeron Fusion contribution, but it does not exhaust the full correction, 
and, in future steps, more complex topologies of the interacting Pomerons have to be considered. 
Moreover, we have approximated the BK evolution kernel to be local in the impact parameter plane. 
In momentum space, this amounts to restricting all Pomerons to carry zero momentum transfer: such 
an approximation can be justified for the scattering on a sufficiently large nucleus; in the case of the 
proton, however, it looks somewhat questionable. Finally, we neglected the contribution from low 
mass excitations of the target proton in the intermediate state, which somewhat reduces the estimated 
magnitude of the rescattering correction. 

Returning to the selection of diagrams, in our previous choice of contributions one recognizes 
a strong asymmetry between projectile and target protons, and it is obvious that a whole class of 
contributions is missing. In physical terms, the Pomerons that give rise to the BK fan diagrams 
propagate in the colour field of the projectile proton and should be screened, too (examples are shown 
in Fig. I^J). It is expected that these new corrections will continue to come with alternating signs, i.e. 
the first screening correction of the screening Pomeron enters the amplitude with the positive sign 
etc. Subsequent iterations of this kind lead to a whole tower of corrections. For the total cross-section 
of nucleus-nucleus collisions an effective theory of interacting BFKL Pomerons |41 | I42[ 14*3]. beyond 
the fan diagram approximation, has been formulated. Imposing a semiclassical approximation that 
resums Pomeron diagrams exemplified in Fig.^, and restricting all Pomerons to have zero momentum 
transfer, numerical solutions have been found. 

All these corrections can be derived from an effective action for interacting Pomerons of Refs. |411 
1421 I4.3| . Thus, a natural continuation of the present study could be an analysis of the hard rescat- 
tering in this framework. At present an evaluation of the interacting Pomeron diagrams of the type 
exemplified in Fig. EH, analogous to the treatment of nucleus-nucleus scattering, might be realistic, 
whereas the diagrams with closed Pomeron loops (see Fig. 02) ) are, to our understanding, still beyond 
reach. Nevertheless, a careful treatment of the effective interacting BFKL Pomeron field theory in 
the semi-classical limit should be capable to provide a more reliable estimate of the absorption of the 
hard gluonic ladders than the one performed in the present paper. 

An important point in our analysis is the role of the saturation scale. Up to now, the analysis 
of diffractive Higgs production had been based upon a rather strict separation between hard and 



23 



a) b) 




Figure 9: Examples of diagrams of interacting BFKL Pomeron field theory that contribute to the 
exclusive Higgs boson production: a) a diagram, without closed Pomeron loops, belonging to the class 
of diagrams that were studied in Refs. \41[ \4%\ \4^j i n the case of nucleus-nucleus scattering; b) a 
diagram with a closed Pomeron loop. 

soft contributions: the hard Two Pomeron Fusion amplitude and the soft rescattering, encoded in 
the survival probability factors. In the first part of our hard rescattering analysis, where a single 
BFKL Pomeron had been considered, we found a strong sensitivity to the infrared scale ko which, at 
first sight, seemed to invalidate the BFKL Ansatz for the rescattering. However, the existence of a 
saturation scale, which for x values being as small as x ~ 10 -6 lies above the infrared region, removes 
this dominance of the small-momentum region and hence supports our perturbative approach. In 
other words, for the production of a Higgs boson which, at the production vertex, provides the large 
momentum scale of Mh, the perturbative nature of the scattering amplitude seems not be limited 
to the production vertex but seems to 'spread out' further away into the rescattering part of the 
amplitude. 

8 Summary 

In this paper we have investigated hard rescattering corrections to the exclusive Higgs boson production 
at the LHC. In a first step, the absorptive correction was given by the exchange of a single BFKL 
ladder which couples to the hard gluonic ladder through the triple Pomeron vertex. This amplitude 
has turned out to be strongly sensitive to the infrared region. In our main part, we then have replaced 
the single BFKL ladder by a Balitsky-Kovchegov gluon density. Due to emergence of a saturation 
scale in the perturbative domain, the infrared behavior stabilizes. We have found that the dominant 
contribution to the hard rescattering correction comes high mass diffractive excitations of the the 
projectile proton, indicating a clear separation of this contribution from that of low mass intermediate 
states of the standard soft rescattering. 
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The correction was found to be negative, with the absolute value being close or even exceeding the 
leading Two Pomeron Fusion contribution. This implies that the effect is significant. However, a more 
precise estimate of the cross sections is not possible until yet more screening corrections are evaluated 
and certain approximation can be lifted. Nevertheless, our results seem to question the simple picture 
of the pp — > pHp process, being described by two hard gluonic ladders that fuse into the Higgs 
boson and leaving space only to a soft rescattering of the protons. We have also discussed possibilities 
of improving, within the framework of the interacting BFKL Pomeron field theory, the description of 
the exclusive Higgs production cross section. 

In order to arrive at a numerical estimate of the exclusive Higgs boson production cross section at 
the LHC, we have calculated the soft gap survival probability for the scattering amplitude containing 
both the Two Pomeron Fusion term and the rescattering correction. For this we have used the 
dependence of the amplitude on the transverse position; due to the different shapes of the Two 
Pomeron Fusion term and of the correction, the gap survival factors were found to depend on the 
strength of the triple Pomeron coupling. 

In our final numerical estimate of the exclusive Higgs boson production cross section we found 
that the cross section may be smaller, and its theoretical uncertainty may be larger than previously 
expected. 

The method developed in this paper may be also applied to other exclusive diffractive production 
processes in pp or pp collisions at the Tevatron or at the LHC. As examples, let us mention the exclusive 
production of di-jets or the exclusive production of heavy scalar mesons in the central region. We also 
view the results of this paper as a step towards developing a better understanding of the effective field 
theory of interacting BFKL Pomerons and of its applications to high energy pp and pp collisions. 
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